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BOUNDARY CONDITIONS FOR HEAT TRANSFER 
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A method for calculating the boundary conditions for heat transfer through glassmaking furnace brickwork is 
presented. The average heat-transfer coefficients and heat fluxes for the principal structural elements of a fur¬ 
nace are determined. A database for mathematical modeling the thermal operation of furnaces is formulated. 
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The coefficient of heat transfer K, W/(m 2 • K), which es¬ 
tablishes a unique relation between the resulting heat flux 
density q, W/m 2 , and the interior temperature f jnl , °C, is a 
given in external heat transfer calculations performed by the 
resolvent zonal method for the surface zones of the brick¬ 
work in the working space of a furnace [1]. 

Mixed boundary conditions are used to determine the 
heat-transfer coefficient. Boundary conditions of the first 
kind are set by the temperature f int on the interior surface of 
the brickwork. Boundary conditions of the third kind are set 
on the exterior surface of the brickwork by the total heat- 
emission coefficient a ext , W/(m 2 • K), and by the tempera¬ 
ture of the ambient medium / amb = 40°C. 

Starting from the formulation of the boundary conditions 
the equation of heat transfer through multilayer brickwork 
has the form 

<7 = ^(4,t-4mb)> (D 

where 


n 1 

z«,+— 

i'=l ® ext 

It is the thermal resistance of the z'th layer of the brickwork, 
(m 2 • K)/W, and n is the number layers in the brickwork. 

In turn 
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where S ; is the thickness of the z'th layer of the brickwork, m, 
and X. (/) is the temperature dependence of the thermal con¬ 
ductivity of the ith layer of the brickwork, W/(m • K). 

For natural convection on the exterior surface of the 
brickwork the coefficient a ext of total heat emission is calcu¬ 
lated from the equation 
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( 4 ) 


where the coefficient & surf takes account of the position of the 
surface and direction of heat emission in space: we take 
k surf = 2.4 for the side walls of the working space and the 
melting tank of the glass furnace (vertical surfaces) and 
k sw[ = 3.3 and 1.6, respectively, for the horizontal surfaces 
with heat emission directed upward and downward (the fur¬ 
nace roof and tank bottom); s surf is the emissivity of the ex¬ 
ternal surface of the brickwork: s surf =0.8 for brick brick¬ 
work and s surf =0.8 and 0.2, respectively, for steel sheet 
and surfaces covered with bright aluminum paint; t ext is the 
exterior temperature of the cladding, °C; and, ct 0 = 
5.67 x 10 ~ 8 W/(m 2 • K 4 ) is the Stefan - Boltzmann con¬ 
stant. 

The gird of the side walls of the melting tank is forced air 
cooled. For f ext = 100 - 400°C the overall heat-emission co¬ 
efficient is given by the expression 

a ext = (9.5 + 0.07f ext )(l+2w), (5) 

where w is the velocity of the cooling air, m/sec. We take 
w = 50 m/sec [2], 
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Fig. 1. Fragment of the transverse section of a glassmaking furnace: 
1 ) DSU brick — 0.4 m; 2) STELLAMUR GLS mastic (0 - 0.3 mm 
grain) — 0.01; 3) LEGRAL 55/05 brick — 0.114 m; 4 ) STELLIT 
GH mixture (0 - 1 mm grain) — 0.030 m; 5 ) MKRV-200 mullite- 
silica felted fabric — 0.060 m; 6) KPD-400 foamed diatomite brick 

— 0.114 m; 7) STELLAMUR GLS mastic — 0.005 m; 8) DSU — 
0.5 m; 9) Z-64 — 0.4 m; 10) RESIMUR ZM 322 mastic (0- 
0.5 mm grain) — 0.005 m; 11 ) MD-1450 — 0.114 m; 12) AZS-33 
block with normal contraction cavity — 0.2 m; 13 ) MD-1650 brick 

— 0.114m; 14) MD-1450 brick — 0.181m; 15) KPD-400 
foamed diatomite plate — 0.067 m; 16 ) AZS-33, NC — 0.5 m; 
17) MD-1450 — 0.067 m; 18) AZS-37 block with reduced contrac¬ 
tion cavity — 0.25 m; 19) MD-1450 — 0.183 m; 20) KPD-500 
foamed diatomite brick — 0.116 m; 21) molten glass — 
1.17/1.60 m (melting/fining zone, respectively); 22 ) AZS-33 plate 
with no contraction cavity- 0.12 m; 23 ) RESISTIT ZM 260 mixture 
(0-3 mm grain) — 0.03 m; 24) Z-64 — 0.065 m; 25 ) BN-60 plate 

— 0.1 m; 26) BN-40 block — 0.2 m; 27) MD-1450 — 0.13 m; 
28 ) KPD-500 — 0.130/0.065 m (melting/fining zone, respectively). 


Taking account of the temperature dependence of the 
thermal conductivity of the brickwork materials and the 
overall heat-emission coefficient, the problem of determin¬ 
ing the heat flux into the surrounding medium is solved by 
the method of successive approximations (method of itera¬ 
tions). The iterative process is stopped when 


density from the exterior surface of the brickwork into the 
surrounding medium; q thT is the density of the heat flux pass¬ 
ing through a wall with thermal conductivity 

_ ^int — ^ext 
<7thr “ > 


where 


s err = 1.0x10 


IS 

-3 




the prescribed computational error, 


An accurate calculation of heat transfer through a multi¬ 
layer wall presupposes a validated temperature given on the 
interior surface of the brickwork. As a rule, the heating tem¬ 
perature of surfaces is taken from experimental data, which 
in the first place are too sparse and in the second place can 
have a very arbitrary relation to the furnace being designed. 
In this connection it is desirable to use the results of numeri¬ 
cal modeling of external heat transfer with optimal plume 
length [3], 

The melting part of the furnace is divided into separate 
structural elements. The roof and longitudinal and end walls 
are singled out in the working space. The melting tank is cal¬ 
culated in greater detail. This pertains to the side walls, 
which are divided along the height into three computational 
sections, and the bottom of the tank [4], The tank is divided 
along its length into a melting zone (70% of the tank length) 
and a fining zone. The forced air cooled surfaces are also sin¬ 
gled out. The structure of the brickwork (see Fig. 1) and the 
characteristics of the refractory and heat-insulating materials 
are prescribed (Table 1) for each surface. 

The nonuniformity of the temperature field on the sur¬ 
faces of the working space and the identical structure of the 
brickwork in the structural elements which have been singled 
out give a basis for using in the calculation of thermal con¬ 
ductivity the average interior temperature ( int . In application 
to the structure of a melting tank with a horseshoe flame the 
values of ( int can be determined from two-dimensional equa¬ 
tions approximating the temperature distribution on the sur¬ 
faces of the working space of the furnace [3]: 

- crown in the working space 

t cr Ofi ? Ti ) = 1 399.65 + 51 L682*i -279.713 y 1 - 

327.564xf + 52049*! y, + 82726 yj 2 ; (6) 

- right-hand longitudinal wall (heated side of the fur¬ 
nace) 


long, rh 


(x^zj) = 1313.145+ 761654*! +490.176zi - 
542979*i 2 -275.787*iz x -1646.567zi 2 ; (7) 


0 ext 0 thr 

q 


<e„ 


- left-hand longitudinal wall (outlet side of the furnace) 


where q is the heat flux density determined from the rela¬ 
tion (1); q ext = a ext (t ext - / amb ) is the density of the heat flux 


hong.ih (*i , z i) = 1278.097+422825*1 +408.297zj - 

208.09lx 2 - 128.85*i zi -1567.367z 2 ; (8) 
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TABLE 1. Thermal Conductivity and Working Temperature of Refractory and Heat-Insulating Articles 


Article brand 

Working tempe¬ 
rature, °C 

Thermal conductivity, W/(m • K) 


DSU, GOST 3910-75 

1600 

0.93 + 0.707 x 10“ 3 f 


STELLAVUR GLS (95.5% Si0 2 ) 

1680 

- 


STELLIT GH (85.5% Si0 2 ) 

1400 

0.29-0.1 x 10 _3 t + 0.3 x 10 _6 t 2 


LEGRAL 55/05 (90% Si0 2 , p = 0.59 g/cm 3 ) 

1500 

0.15 + 0.2 x 10 _3 t- 0.04 x 10 _6 t 2 


MKRV-200, GOST 23619-79 

1150 

0.04 + 0.1 x 10 _3 t+ 0.08 x 10“ V 


KPD-400, TU 5764-002-25310144-99 

950 

0.095-0.757 x 10“ 3 t+0.25x 10“ 6 t 2 


Z-64 (64% Zr0 2 ) 

1630 

2.7-0.01 x 10“ 3 t +0.077 x 10' 6 t 2 


MD-1450 (52% A1 2 0 3 , p = 0.8 g/cm 3 ) 

1430 

0.23 + 0.1 x 10“ 3 t 


MD-1650 (73% A1 2 0 3 , p = 1.0 g/cm 3 ) 

1430 

0.384 + 0.07 x 10 _3 f 


AZS-33 (33% Zr0 2 ) 

1700 

8.84-11.9 x 10~ 3 t + 7.0 x 10 “V 


AZS-37 (37% Zr0 2 ) 

1700 

6.0-5.63 x 10 -3 t + 3.86 x lO -6 ^ 


BN-40 (43% A1 2 0 3 ) 

1350 

1.02 + 0.62 x 10~ 3 f- 0.3 x 10 ~ 6 t 2 


BN-60 (62% A1 2 0 3 ) 

1550 

1.379 + 0.32 x 10“ 3 f 


RESISTIT ZM 260 (18% Zr0 2 ) 

1650 

2.5- 1.8 x 10 _3 t + 0.8 x 10“ 6 t 2 


RESIMUR ZM 362 (31% Zr0 2 ) 

1650 

5.0-5.1 x 10 3 f + 2.32 x 10“¥ 



- end wall (flow-through) 

f e .fi ( Vi ) = 1575.210-14.30lj! -350.668zi + 

24.909^ -36L626y 1 z I -398.718zj 2 , (9) 

where x x , y v Zj are, respectively, x/L { , y/L f and z/T f ; L f is 
the length of the working space of the furnace, m. 

Averaging Eqs. (6) - (8) gives: t a = 1477.9°C, f long rh = 
1526.7°C, t long lh = 1434.0°C and t c fl = 1526.7°C, respec¬ 
tively. For the inlet end wall we set t e m = 1355.2°C [3]. 

Working space of the furnace. We shall now examine 
the computational results for the parameters of heat transfer 
through the brickwork of the working space of a furnace. 

Crown (Table 2). Modern heat insulation is installed on 
the main part of the surface of the crown. The low value of 
the heat-transfer coefficient 0.49 W/(m 2 • K), due to the high 
thermal resistance of the brickwork 1.99 (m 2 ■ K)/W, makes 
it possible to decrease the heat losses into the ambient me¬ 
dium from 4786.4 to 705.4 W/m 2 , i.e., 6.8-fold. 

The expansion joints between crown sections and the in¬ 
terlocking row of the brickwork are left uninsulated because 
of safety considerations. The area of the uninsulated surface 
is 7.44% of the crown area. Even though the heat losses 
through the uninsulated surface of the crown are high, the av¬ 
erage coefficient of heat transfer K is characterized by a low 
value 0.70 W/(m 2 • K), while the average heat flux q does 
not exceed 1009.0 W/m 2 . 

Longitudinal walls of the working space (Table 3). The 
computational results show that the average heat losses 
through the longitudinal walls depend strongly on the con¬ 
struction of the walls (see Fig. 1). About 40% of the area of a 
longitudinal wall falls on sections with uninsulated (I) and 
inadequately insulated refractory brickwork (II, IV). As a re¬ 


sult, if the heat fluxes through the heat-insulated section of 
the heated and outlet sides of the furnace are 968.7 and 
892.5 W/m 2 , respectively, then the average value q for the 
indicated sides reaches 2104.0 and 1916.3 W/m 2 . 

A detailed accounting of the structural features of a lon¬ 
gitudinal wall makes it possible to determine a realistic value 
of the average coefficient of heat transfer. This value is > 2 
times greater than the value of K for an insulated section of 
the wall. 

End walls of the working space (Table 4). An end wall 
consists of a thermally insulated section (III) and a bearing 
baddeyelite-corundum tooth (IV). The structure of the ther¬ 
mal insulation of both sections is similar to that of a longitu¬ 
dinal wall (see Fig. 1). 

The areas of the heat-insulated surface of the flow¬ 
through and inlet walls are 88.6 and 77.7%, respectively. In 
consequence, the average heat-transfer coefficient and heat 
flux are much lower than for a longitudinal wall. Neverthe- 


TABLE 2. Computational Results for Heat Transfer through the 
Furnace Crown 


Computational 

parameters 

Insulated section 
of the crown 

Uninsulated section 
of the crown 

7 =7 °r 

hnt *cr> ^ 

1477.9 

1477.9 

t °C 

*ext» ^ 

83.6 

229.9 

R, (m 2 • K)/W 

1.99 

0.26 

“ext> W/(m 2 ■ K) 

16.17 

25.23 

K , W/(m 2 ■ K) 

0.49 

3.33 

q, W/m 2 

705.4 

4786.4 

K , W/(m 2 • K) 

0.70 

q, W/m 2 

1009.0 
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TABLE 3. Computational Results for Heat Transfer through a Longitudinal Wall of a Furnace (Heated/Outlet Sides) 


Computational parameters 


Section of the wall surface 


1 

ii 

III 

IV 

^int — ^long.rhAint — ^long.lh’ ^ 

1526.7/1434.0 

1526.7/1434.0 

1526.7/1434.0 

1526.7/1434.0 

t °C 

*ext > ^ 

226.1/215.7 

169.2/162.8 

106.9 /102.7 

224.2/211.6 

R, (m 2 • K)/W 

0.32/0.33 

0.58/0.58 

1.46/1.49 

0.33 /0.34 

a ext , W/(m 2 ■ K) 

21.64/20.98 

18.14/17.76 

14.48/14.22 

21.52/20.72 

K, W/(m 2 ■ K) 

2.71/2.64 

1.58/1.56 

0.65/0.64 

2.67/2.55 

q, W/m 2 

4030.2 /3686.6 

2344.2/2181.5 

968.7/892.5 

3966.2/3552.0 

K, W/(m 2 • K) 


1.42/1.37 


q, W/m 2 


2104.0/1916.3 



less, it should be noted that q and K differ considerably from 
the values characterizing the heat-insulated sections of the 
walls. 

Melting tank of the furnace. We shall now examine the 
computational results obtained for the parameters of heat 
transfer through the brickwork of a melting tank filled with 
dark-green molten glass. The temperature dependence of the 
effective thermal conductivity L eff of the molten glass is re¬ 
presented by the equation 

L eff = 1.0908-0.5 x 10' 3 ? + 6 x 10- 6 t 2 . (10) 

Since heat is transferred into the molten glass through its 
exterior surface, the main boundary condition for the calcula¬ 
tion is the distribution of the temperature t mg on this surface. 
In two-dimensional dimensionless coordinates we have [3] 

t mg (x u y x ) = 1132904 + 93L278VJ +67.537^ - 

675.024* f + 4.71&c 1 y 1 -130.010^, (11) 

where x v y x are x/L mt and y/L m , respectively, and L mt is the 
length of the melting tank, m. 


TABLE 4. Computational Results for Heat Transfer through the 
End Walls of a Furnace 


Section of wall surface 


Computational 

parameters 

Inlet wall 

Flow-through wall 

hi 

IV 

hi 

IV 

^int — 4. ft Aint — 4. in’ ^ 

1355.2 

1355.2 

1525.9 

1525.9 

'ext-°C 

99.2 

201.2 

106.8 

224.1 

R, (m 2 ■ K)/W 

1.51 

0.35 

1.46 

0.33 

a ext’ W /(m 2 'K) 

14.01 

20.07 

14.47 

21.51 

K, W/(m 2 • K) 

0.63 

2.46 

0.65 

2.67 

q, W/m 2 

829.7 

3233.6 

968.0 

3962.4 

K, W/(m 2 • K) 

1.04 

0.88 

q, W/m 2 

1364.8 

1309.1 


Let us transform Eq. (11) into the one-dimensional fonn 

t^"(x 1 ) = 1236.3-30.112*! + \21A.6xl - 

613.08.vj 5 -539.06V! 4 , (12) 

where i mg is the average temperature of the molten glass sur¬ 
face over the width of the tank, °C. 

We obtain from Eq. (12) the average values of the sur¬ 
face temperature of the molten glass for the melting and fin¬ 
ing zones: f mg = 1369.1 and 1415.3°C, respectively; we shall 
take these values as the boundary conditions for calculating 
heat transfer through the melting tank. The molten glass is 
treated as the first layer of a multilayer wall. 

Tank bottom (Table 5). The computational results show 
that the proposed structure of the brickwork at the bottom of 
the tank gives it a high thermal resistance R =1.88 and 
1.43 (m 2 • K)/W and, in consequence, a low value of the 
heat-transfer coefficient W=0.51 and 0.66 W/(m 2 • K) for 
the melting and fining zones, respectively. The effectiveness 
of the thennal insulation of the bottom brickwork gives low 
heat losses into the ambient medium ( q = 679.5 and 
914.1 W/m 2 ), as a result of which the temperature of the 
melt layer at the bottom is high ? bot = 1297.4 and 1286.6°C, 
corresponding to the conditions for high-capacity glass¬ 
making. 


TABLE 5. Computational Results for Heat Transfer through the 
Melting Tank of the Furnace 


Computational 

parameters 

Melting zone 

Fining zone 

7. =t °C 

hnt mg’ 

1369.1 

1415.3 

Wc 

1297.4 

1286.6 

Lxt ’ ° c 

98.0 

113.0 

R, (m 2 • K)/W 

1.88 

1.43 

«ext> W/(m 2 ■ K) 

11.72 

12.51 

K, W/(m 2 • K) 

0.51 

0.66 

q, W/m 2 

679.5 

914.1 
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TABLE 6. Computational Results for Heat Transfer trough the Side Walls of the Melting Tank of the Furnace 




Melting zone 



Fining zone 


Computational 

parameters 


Wall section 



Wall section 



V 

VI 

VII 

V 

VI 

VII 

f. °p 

1369.1 

1347.5 

1316.6 

1415.3 

1368.4 

1311.9 

t °C 

*ext» ^ 

134.0 

99.6 

98.0 

136.6 

100.6 

97.8 

R, (m 2 • K)/W 

0.06 

1.50 

1.51 

0.06 

1.49 

1.51 

a ext , W/(m 2 • K) 

207.68 

14.03 

13.91 

209.66 

14.09 

13.92 

K, W/(m 2 ■ K) 

14.68 

0.64 

0.63 

14.72 

0.64 

0.63 

q, W/m 2 

19512.3 

836.0 

808.0 

20235.1 

854.0 

804.1 

K, W/(m 2 • K) 


2.91 



2.56 


q, W/m 2 


3848.5 



3478.8 



Side walls of the melting tank (Table 6). We shall take a 
linear temperature distribution over the thickness of the layer 
of glass. Then it is possible to determine the average heating 
temperature of the surfaces of the side walls of the tank at 
sections with discrete and continuous insulation [4]. For the 
cooled sections of the walls in the melting and fining zones 
we take t mt = t m% = 1369.1 and 1415.3°C, respectively. 

It follows from the data in Table 6 that even though the 
heat insulation of the walls is effective the heat losses 
through the cooled section have a dominating influence on 
the average heat-transfer coefficient and heat flux into the 
ambient medium. Since cooling intensity is dictated by the 
technological aspects of the work of the refractories [2], the 
heat losses through the gird of the tank can be seen as justi¬ 
fied. Hence the values K = 2.91 and 2.56 W/(m 2 • K) (melt¬ 
ing and fining zones, respectively) reflect the real conditions 
of heat transfer through the side walls of the melting tank. 

In closing, we note that the values presented in this arti¬ 
cle for the average heat-transfer coefficients can be regarded 
as representative initial data for mathematical modeling of 


the thermal work of high-capacity glassmaking furnaces. The 
computational results for the heat losses through the brick¬ 
work of the working space and melting tank can be used for 
constructing the heat balance of glassmaking furnaces with 
the energy-conserving structure of modern furnace enclo¬ 
sures. 

REFERENCES 

1. V. Ya. Dzyuzer, V. S. Shvydkii and E. B. Sadykov, “Numerical 
model of external heat transfer in a gas-electric glassmaking fur¬ 
nace,” Steklo Keram., No. 2, 8 - 12 (2012). 

2. V. Ya. Dzyuzer, “Cooling of the gird of the melting tank of a 
glassmaking furnace,” Steklo Keram., No. 8, 13 - 16 (2009). 

3. V. Ya. Dzyuzer and V. S. Shvydkii, Designing of Energy-Efficient 
Glassmaking Furnaces [in Russian], Tcplotekhnik, Moscow 
(2009). 

4. V. Ya. Dzyuzer, “Designing the refractory brickwork of the melt¬ 
ing tank of high-capacity glassmaking furnaces,” Ogenup. Tekh. 
Keram., No. 8, 14 - 22 (2007). 








